Kondo proximity effect: How does a metal penetrate into a Mott insulator? 
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We consider a heterostructure of a metal and a paramagnetic Mott insulator using an adaptation 
of dynamical mean field theory to describe inhomogeneous systems. The metal can penetrate into 
the insulator via the Kondo effect. We investigate the scaling properties of the metal-insulator 
interface close to the critical point of the Mott insulator. At criticality, the quasiparticle weight 
decays as 1/x 2 with distance x from the metal within our mean field theory. Our numerical results 
(using the numerical renormalization group as an impurity solver) show that the prefactor of this 
power law is extremely small. 

PACS numbers: 73.20.-r, 71.27.+a, 71.30.+h 



Introduction: In the last few years, an enormous 
amount of interest has arisen in heterostructures fabri- 
cated out of strongly correlated materials. Driven by the 
prospect of new effects and devices based on correlated 
electron compounds, a wide range of systems has been 
studied experimentally and theoretically. For example, 
the interface of two Mott insulators can show metallic be- 
havior 0, [3, H, IH HI or can even become superconducting 
[(| . The conducting layers can, for example, arise from 
interface charges induced by the Coulomb interactions. 
As shown by Thiol et al. |5| such interfaces can be ma- 
nipulated by gate voltages thereby opening the prospect 
for interesting novel devices. 

The physics of such inhomogeneous systems can also 
play a role for the properties of bulk materials where 
competing phases lead to the formation of domain walls. 
For example, it has been argued in Ref. [7| that the con- 
ductivity close to the endpoint of the first-order Mott 
transition in certain organic salts of the k-ET family is 
dominated by domain wall effects. Also for cold atoms, 
the trapping potential naturally makes the experimental 
systems inhomogeneous which often leads to coexisting 
phases and corresponding phase boundaries Q . 

In this paper we will investigate the interface of a metal 
and a Mott insulator. How does a metal penetrate into 
a Mott insulator? The main difference between a Mott 
insulator and an ordinary band insulator is the presence 
of magnetic degrees of freedom arising from the localized 
spins. While the large charge gap, of the order of the 
local Coulomb repulsion U, prohibits tunneling of elec- 
trons into a Mott insulator, the resonant spin flip scatter- 
ing opens a new channel for tunneling via the well-known 
Kondo effect and allows metallic behavior to be induced 
within the Mott insulator. Due to this 'Kondo proxim- 
ity effect' an insulating layer adjacent to the metal will 
also become metallic. In this manner the metal 'eats' it- 
self layer by layer into the Mott insulator if not stopped 
either by magnetism or thermal fluctuations. Here, we 
study this physics within the simplest setup consisting 
of a particle-hole symmetric Hubbard model (see below) 



where the local interaction U jumps across the interface 
from U = C/icft <C U c to a value U = [/right close to the 
critical coupling U c of the Mott transition. For such a 
model the charge is always homogeneous and no compli- 
cations due to charge-reconstruction or charged surface 
layers arises. Furthermore, we only consider a paramag- 
netic Mott insulator and comment on the role of mag- 
netism only in the conclusions. 

For one-dimensional systems, powerful numerical 
methods like DMRG are available to study inhomoge- 
neous strongly interacting models jsj. For three dimen- 
sional systems, however, further approximations are nec- 
essary. Here, the method of choice to study the Mott 
transition is the so-called dynamical mean field theory 
(DMFT) 0, El. Within DMFT, the only approxi- 
mation is to neglect non-local contributions to the self- 
energy. This app roximation can be used both for homo- 
geneous [13, U| and inhomogeneous [H, [H, 0, 0, EH 
problems. In the case of a heterostructure, each layer is 
effectively mapped to a single-impurity Anderson model. 
These are coupled by a self-consistency condition, see be- 
low. 

Potthoff and Nolting used this spatially resolved 
DMFT to study the Mott transition at surfaces using 
the semi-infinite Hubbard model [l6j |. In this context, 
also the question was studied how a metallic surface influ- 
ences the insulating bulk. Spatially resolved DMFT has 
also been applied [l[ [H, 13| to investigate heterostructures 
of Mott insulators and band insulators. 

A main problem of DMFT is the need for a reliable and 
efficient method to solve the effective impurity problem. 
Previous applications of DMFT to inhomogeneous sys- 
tems used impurity solvers based on exact diagonaliza- 
tion of small systems [3] , a linearized version of DMFT 
close to the critical point [16], a two-site approxima- 
tion [J [j| or slave-boson mean-field theory [ijj , implying 
severe further approximations, or started from simpler 
models such as the Falicov- Kimball model [3]. Only 
recently [H, HH, the numerical renormalization group 
(NRG) method [H HJ] was implemented as an impu- 
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rity solver to study the Mott transition of trapped atoms 
in an optical lattice. We will also use this approach here 
as the NRG appears to be the only method presently 
available which can quantitatively resolve quasiparticle 
weights as small as 10~ 3 which are needed to describe 
the physics close to the Mott transition. 

After introducing the model and our method (DMFT 
for inhomogeneous systems+NRG), wc will first investi- 
gate the heterostructure at finite T and for £/ r i g ht = U c . 
We will then analyze the T — scaling properties of the 
interface region at C/ r i g ht < U c and t7 r ight > U c using both 
DMFT+NRG and a Ginzburg-Landau type analysis. 

Model and Method: To investigate the junction of the 
metal and the Mott insulator, we will consider the half- 
filled Hubbard model 
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FIG. 1: (Color online) Layer dependence of local spectral 
function close to the interface x = for C/ r i g ht = U c and 
T = 1.14 x 10" 5 D. Inset: A(u) near uj = 0. 



on a three-dimensional cubic lattice with the half band 
width D — 6i. While we will consider a uniform hopping 
t, we choose U — U\ c a = D for x < describing a metal 
with a sizable quasiparticle weight Z mcta i — 0.62. For 
sites with x > 1, we use an interaction ?7j = f7 r i gri t ~ U c 
close to the the critical value, U c « (2. 79 ±0.01)1? which 
separates the metallic from the insulating phase in the 
bulk. Note that the system is translationally invariant in 
the yz directions as Ui is constant within each yz layer. 

The DMFT algorithm for this heterostructure is al- 
most identical to the standard one [lo| . As all sites 
within a single yz layer arc equivalent, it is sufficient 
to solve only one effective Anderson impurity problem 
(using NRG [2(| [2l[) for each yz layer to obtain an x 
dependent self energy T, x (oj). From this one obtains the 
lattice Greens function 

G lat (e ki ,^) = (w-ek x -t xx , -.WSxM)- 1 (2) 

written as a matrix in the x coordinates where t xx i = t 
for x' = x ± 1 and otherwise and e ki is the dispersion 
within each layer. From G lat one determines the local 
Greens function which is used (loj to derive a new effec- 
tive Anderson impurity model with the Greens function 
Gi m p (a;) for each layer using the self-consistency equation 



G7 P (w) = / de k± N 2d (e k± ) G lat (e k ^ 



(3) 



where N 2 d{t) = J2k ± — £ kj_) is the two-dimensional 
density of states of the yz layers. To avoid numeri- 
cal difficulties associated with the logarithmic divergence 
of 2d cubic density of states, we use A^2d(e) = l/(8£) 
for |e| < At. This does not affect any universal prop- 
erties discussed below and leads only to small changes 
in U c (and other numerical prefactors) of about 10%. 



From Eq. ^ it seems that one has to invert a ma- 
trix for each value of u> and ek^ • Fortunately, this com- 
putationally expensive step can be simplified by diago- 
nalizing G lat (0, w)^ 1 using the orthogonal matrix 0(u>), 
O T (uj)G lat (0,u)- 1 d(w) = M(w), so that G lat (e ki ,w) = 
6(u)(M(u) - e k± )- 1 6 T (uj). We use 20 metallic layers 
with U = D and 40 layers with U ~ U c which is suffi- 
ciently large to avoid any finite size effects. 

Results: Fig. Q]shows the layer dependence of the spec- 
tral function for t7 r ; g ht — U c for x > and for a low 
temperature T — 1.14 x 10~ 5 Z3. The metallic side is 
only weakly affected by the presence of the insulator as 
within our paramagnetic, particle-hole symmetric model 
no Friedel oscillations occur. All layers with U — U c show 
pronounced Hubbard bands. The width of the sharp 
quasiparticle peak, which describes the penetration of 
the metal into the quantum critical Mott state, decays 
rapidly. The quasiparticle peak collapses completely (up 
to an exponentially small feature) from the 5th layer on, 
when the Kondo temperature of the corresponding im- 
purity model becomes much smaller than T. 

For a quantitative analysis of how the metal penetrates 
into the Mott insulator we investigate a heterostructure 
consisting of a 'good metal', U\ it = D and a 'bad metal', 
bright < C/ c at T = 0. For T = the quasiparticle weight 
Z x of layer x is well defined and can be obtained from 
Z x = (1 — <9 w Re£a;(ti;)) -1 . Fig. [2] shows the quasipar- 
ticle weight Z as a function of the distance x from the 
interface. Upon increasing C/ r i g ht the quasiparticle weight 
deep in the bad metal decreases linearly with U c — t7 r j g ht 
[Ioj |. Close to the critical point one expects scaling be- 
havior and indeed we observe in Fig. [3] 
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with v = 1/2 where f[u] is an universal scaling function 
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FIG. 2: (Color online) Quasiparticle weight Z for an inho- 
mogeneous layered system describing a 'good' metal (f = 
Cleft = D for x < 1) in contact with a 'bad' metal (fright < 
U c = (2.79 ± 0.01)13 for x > 1). Inset: Quasiparticle weight 
for a heterostructure of a good metal and a Mott insulator 
(Kight > Uc)- In both cases, Z drops as 1/x 2 for x < see 
Fig. [3] For extremely small quasiparticle weights Z < 10~ 3 
the numerical results eventually become unreliable. 

with f[0] = 1 and f[u -> oo] » (0.150 ± 0.005)u 2 for 
bright U c . The observation that DMFT is character- 
ized by the usual mean-field exponent v — 1/2 and the 
1/x 2 decay of the correlation function in the quantum 
critical regime is one of the main results of this paper. 
Defining the correlation length £ by f[u] = 2, we obtain 

f U e \ 1/2 

\U C — U 'right / 

In the Mott insulating phase at T = 0, fright ^ U c , the 
quasiparticle weight drops rapidly, see inset of Fig. [2] In 
this regime, the numerical calculations become difficult 
and are unreliable for Z < 10 ~ 3 . Nevertheless, in the 
regime where the numerical results are accurate, they 
confirm the scaling ansatz Eq. ([J), as can be seen in detail 
in Fig. [3J In the quantum-critical regime, i.e. for x < £, 
the quasiparticle weight decays as 1/x 2 with the same 
prefactor as in Eq. (j4|). For x > £, however, Z drops 
exponentially but remains always finite. 

This picture is further corroborated by an analysis in 
the spirit of a Ginzburg-Landau mean-field treatment as 
in [3, HI)]. The basic idea is that close to U c , the 
physics is mainly determined by the quasiparticle peak 
which can be characterized by a single number, the quasi- 
particle weight Z x . We therefore approximate in Eq. ([2]) 
E x (w) ss lu — oj / Z x . The resulting local spectral function 
— - ImGjj(w) has a peak with a finite width. From this 
peak one has to determine a single number describing the 
effective impurity model and from this a new value of Z x 
using a Ginzburg-Landau expansion around the critical 
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FIG. 3: (Color online) Scaling plot 24] of the T = quasi- 
particle weight close to the quantum critical point both for 
fright < U c (upper curves) and [/right > U c (lower curves) 
confirming v = 1/2 in Eq. [4] In the quantum-critical regime 
(a; < £)> z w 0.008/a; 2 (solid line). For fright < U c , Z 
saturates at a finite value proportional to U c — fright for 
x > ^ while it drops exponentially for f right > U c , see Eq. [8] 
Dashed, dot-dashed lines: scaling curves obtained from a 
Ginzburg-Landau type analysis using Eqs. © and {?). The 
results for ZU C / (fright — U c ) < 0.05 are not numerically reli- 
able. 

point. 

K = ^(Z X _ 1 + ^Z X + Z X+1 ) (6) 

Z x = Z' x - a - - Uc Z' x - (3Z' X 2 (7) 

For the first step, Eq. ([6]), we have used the procedure 
described by Bulla, Potthoff and Nolting [H, HI with 
/ N2d{t)de = 16t 2 /3. For the homogeneous system with 
U <U C one obtains Z = % V \J JJ ■ Analyzing the asymp- 
totic solutions of Eqs. © and ([7]) for the interface with 
fright = U c we find lim^^oo Z x — 9/(ll(3x 2 ). The model 
(|6l7p reproduces the critical exponents of our DMFT cal- 
culation (|4]) . Fitting the asymptotic formulas to the NRG 
results we obtain a = 15.7 ± 5 and (3 = 102 ± 30, reflect- 
ing the small prefactor in Eq. ((4}. While Potthoff and 
Nolting [lfl have analyzed different critical exponents, 
their results are qualitatively fully consistent with ours. 

In Fig. [3l we show that the Ginzburg-Landau analysis 
reproduces the DMFT scaling curves quite well as the 
parameters in I|6l7p have been determined to fit the pref- 
actor of the 1/x 2 law and the large x limit for L^ght < U c . 
The large deviations for fright > U c arise in the afore- 
mentioned regime where the numerical calculations are 
no longer reliable. There is also a smaller deviation for 
fright < U c which could be an indication that the model 
(|6I7|) does not reproduce DMFT even in the scaling limit 




4 



|25l |. From the model (|6l7p one can also extract the 
asymptotic behavior for large x > £ 

0.8 U r 
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0.09 



for 



"right — U c 



(8) 
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Conclusions: In this paper, we have studied how a 
metallic state penetrates into a paramagnetic Mott in- 
sulator (or a bad metal). Using a scaling analysis close 
to the quantum critical point we have determined within 
dynamical mean field theory the critical exponents and 
the asymptotic behavior of the quasiparticle weight close 
to and far away from the interface. 

The main physical mechanism governing the interface 
of a metal and a Mott insulator is the Kondo effect: the 
localized spins of the Mott insulator are screened when 
they are brought into contact with the metal and be- 
come therefore part of the metal. However, our numeri- 
cal results show that this mechanism is not very effective 
as can be seen from the numerical prefactor in Eq. ([4]). 
Even for J7 r i g ht = U c , the quasiparticle weight is only of 
size 0.008/x 2 . There is no small parameter in the model 
which controls this prefactor, which is reminiscent of an- 
other small number characterizing the physics of Mott 
insulators: the critical temperature describing the end- 
point of the first-order Mott transition is much smaller 
than the Mott gap both within DMFT and in systems 
like V2O3 |26j. Also the correlation length is extremely 
short: to obtain in Eq. © a correlation length of 10 lat- 
tice spacings, one has to approach the critical point with 
a precision of 10~ 4 . For all practical purposes our results 
imply that the Mott insulator is de facto impenetrable to 
the metal: Mott insulators are very good insulators and 
the 'Kondo proximity effect' is inefficient. This is consis- 
tent with our previous study of trapped fermionic atoms 
in an optical lattice Q, where a metallic phase barely 
penetrates into a coexisting Mott insulator. 

The small quasiparticle weights at T = also imply 
that very small temperatures larger than the local Kondo 
temperature, Tk oc Z, efficiently quench the 'Kondo 
proximity effect', see Fig. [TJ Even more important is 
the effect of magnetism which we have neglected in our 
study. The tiny local Kondo temperatures in the Mott 
insulating phase will typically be much smaller than the 
exchange couplings of the spins, wiping out the Kondo ef- 
fect. The magnetism of the Mott insulating phase will, in 
contrast, penetrate easily into the metal [18j via Friedel 
oscillations of the magnetization. 

For the future, it will be interesting to investigate with 
our methods also models which are not particle-hole sym- 
metric where interface charges and long range Coulomb 
interactions can lead to an electronic reconstruction of 
the interface 
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